From Weakly to Strongly Magnetized Isotropic MHD Turbulence 
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High Reynolds number isotropic magneto-hydro-dynamic turbulence in the presence of large scale 
magnetic fields is investigated as a function of the magnetic field strength. For a variety of flow 
configurations the energy dissipation rate e follows the Kolmogorov scaling e oc even when 

the large scale magnetic field energy is twenty times larger than the kinetic. Further increase of the 
magnetic energy showed a transition to the e oc Ur ma B rms / 1 scaling implying that magnetic shear 
becomes more efficient at this point at cascading the energy than the velocity fluctuations. Strongly 
helical configurations form helicity condensates that deviate from these scalings. Weak turbulence 
scaling was absent from the investigation. Finally, the magnetic energy spectra showed support for 
the Kolmogorov spectrum k~ 5 ^ 3 while kinetic energy spectra are closer to the Iroshnikov-Kraichnan 
spectrum fc~' i//2 . 
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One of the most fundamental questions that can be 
asked about an out-of equilibrium system is the relation 
between the energy injection/dissipation rate e, and the 
amplitude of the fluctuations ug. In hydrodynamic tur- 
bulence such estimates are clear and the desired relation 
comes from the balance between the injection rate and 
the flux of energy to the small scales due to nonlinear 
interactions. Such considerations lead to the turbulent 
scaling e oc uf/£ioi Re cx ug£/v 3> 1. Here v is the 
kinematic viscosity, ug is the amplitude of the velocity 
fluctuations at the forcing scale I = £ F . Constancy of 
the energy flux over all scales I leads to m oc e 1 / 3 ^ 1 / 3 
that results to the Kolmogorov (K41) prediction for the 
energy spectrum E(k) = C K e 2 / 3 k~ 5 / 3 Although the 
proportionality coefficient C K for the turbulent scaling 
is still debatable in the literature @], the scaling is well 
supported by both experiments and numerical simula- 
tions di]. 

The situation becomes more complex when a linear 
wave term is introduced in the dynamical equations. 
Such terms introduce new timescales in the system and 
the validity of the previously derived relations is in ques- 
tion. If the amplitude of the wave frequency u is suffi- 
ciently larger than the amplitude of the fluctuation shear 
ug/£ then the system can be treated within the frame 
work of weak turbulence theory Q . Then depending on 
the wave resonances (three-wave, four-wave etc) the en- 
ergy injection rate is decreased and becomes 
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where n = for strong turbulence and n = 1 for three- 
wave interactions @ ■ As the strength of the wave term is 
increased a transition is expected from the Kolmogorov 
scaling for which the energy dissipation rate scales like 
the fluctuation amplitude to the cubic power to the wave 
turbulent scaling where higher powers are involved. 

Magneto-Hydro-Dynamic (MHD) turbulence is an ex- 
ample of turbulence where eddies and (Alfven) waves co- 
exist [f| , and provides a good approximation to a va- 
riety of industrial and astrophysical flows @, . When 



a flow is coupled to a uniform magnetic field Bq fluc- 
tuations travel parallel or anti-parallel to the magnetic 
field lines with cj oc ±Bq/£\\. The indexes _L, || indi- 
cate the direction perpendicular and parallel to the mag- 
netic field respectively. The nonlinearities allow for three 
wave interactions [l(| and the scaling e = Cu\/Bq£ is 
expected. If isotropy is assumed (£± ~ £n ~ £), con- 
stant energy flux over all scales £ implies ut oc (eBol) 1 ^ 4 
that leads to Iroshnikov-Kraichnan (IK) energy spectrum 
E(k) oc (ei?o) 1/2 fc- 3/2 [HEl. 

However the assumption of isotropy that was assumed 
is not valid in the presence of a strong magnetic field that 
suppresses the cascade in its direction [13J. For strong 
turbulence K41 arguments can be repeated taking into 
account anisotropy that leads to the energy spectrum 
E(k) oc e 2 / 3 fc ± 5 ^ 3 with the so called critical balance re- 
lation ue/£± ~ Bq/£u [Tij . An anisotropic IK energy 

spectrum can also be recovered E(k) oc e 2 / 3 fc ± 2 ^ 3 by as- 
suming that the proportionality coefficient C that ap- 
pears in cq. [1] decreases with scale like C oc i 1 / 4 due 
to scale dependent alignment of the two fields (veloc- 
ity and magnetic) that quenches the non-linearities [l5j |. 
The two proposed spectra have been investigated, in the 
last years by numerical simulations 0, by differ- 
ent groups without however reaching agreement for the 
strong turbulence spectrum. 

Weak turbulence on the other hand valid for ug/£± <C 
Bp/£u can be treated perturbativcly. Detailed calcula- 
tion [HI leads to the prediction for the spectrum E(k) oc 
f{k\\)k^ 2 {eB (s ) 1 / 2 k1 2 . However weak MHD turbulence 
has the particularity that all three-wave resonant inter- 
actions contain one wavenumber with zero projection in 
the direction of the magnetic field and zero phase veloc- 
ity. The set of these wave numbers themselves compose 
a strongly coupled nonlinear system and thus further as- 
sumptions are needed for the validity of weak turbulence. 
In numerical simulations weak MHD turbulence has been 
observed only when the resonant manifold k\\ — is not 
forced [134201 • When it is forced the 2D components 
dominate, the turbulence is strong and can also introduce 
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an inverse cascade as observed in 2-dimcnsional flows [2C§ . 

In nature uniform magnetic fields do not exist. How- 
ever, magnetic fields Bl that vary over large length-scales 
L can be approximated as uniform provided that turbu- 
lent energy remains in much smaller scales. The validity 
of this approximation however is in doubt since small 
scale variations £j_ -c L couple to parallel variations 
£ii ~ Bl£±/u£ that can be as large as L provided that 
Bl is strong enough. Then the assumption of uniformity 
is not valid. Thus, MHD turbulence in the presence of 
large scale (but not uniform) fields needs to be revisited. 

To study statistically isotropic MHD turbulence in the 
presence of large scale magnetic fields we employ high 
resolution direct numerical simulations of the MHD equa- 
tions: 

9 f u + u-Vu = b-Vb-VP + i/V 2 u + F u (2) 
<9fb + u- Vb = b- Vu + r/V 2 b + F b . (3) 

in a triple periodic box of size L = 2ir. Here u is the ve- 
locity field and b the magnetic field. Both fields satisfy 
V • u = V • b = and (u) = (b) = 0, where the angular 
brackets stand for spatial average, v is the viscosity and 
r\ the magnetic diffusivity. F u is an external mechanical 
body force, while Fb = V x £ where £ is an external elec- 
tromotive force. F„ and F;, are both solenoidal functions 
varying randomly in time with time correlation r. F u is 
acting only on wavenumbers with |k| = k u = 2 and is 
non-helical: (F u -V x F u ) = 0. Fb is acting only in the 
largest scale of the system |k| = kf, = 1, and in general 
has non-zero helicity. All the parameters of the runs can 
be found in table I. For the simulations a pseudo-spectral 
code was used [U [22[ on grids of size 512 3 (Runs A#) 
and f024 3 (Runs B#). 

In a typical helical run all quantities grow initially up 
to a point when dissipation rates and kinetic energy reach 
a steady state while the magnetic energy is still increasing 
slowly. During this time the magnetic field is composed 
of a large scale helical component Bl with |k| ~ 1 that 
contains most of the magnetic energy and small scale tur- 
bulent fluctuations bg of amplitude bg ~ u. Thus mag- 
netic energy E M provides a measure of the large scale 
field E M ~ \B\, while kinetic energy E K provides a 
measure of the turbulent fluctuations. The growth of Bl 
depends on the amplitude and the helicity of the mag- 
netic forcing. The time evolution of the large resolution 
runs can be seen in figure [TJ Due to the slow increase of 
the magnetic energy it is possible to perform short time 
averages (over a few turn over times) and thus obtain 
global averaged quantities for various values of magnetic 
energy from the same run. Thus the turbulent scaling of 
the energy dissipation rate can be tested from multiple 
measurements for a wide range of magnetic field strength 
that here it is quantified by fx = E M / E K . The isotropy of 
the system also allows to to perform spherical averages 
and thus improve two point statistics compared to the 
case with uniform magnetic fields. 

Figure [2] shows the energy dissipation rate normal- 
ized by U^ ms k u = (2E K ) 3 / 2 k u as a function of the en- 
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TABLE I. Table with the parameters of all runs. The 
numbers in the first four columns are the input parame- 
ters given by Q = (F u 2 ) 2 /i/ 2 ^ 3 , the Grashof number, M = 
<F b 2 )2/(F u 2 >2 the ratio of magnetic to mechanical forc- 
ing, H = (F b • V x Fb)/(F b 2 )fci, the relative helicity of the 
forcing and St = Tf/r the Strouhal number where Tf = 
l/((F u 2 )fc 2 ) 1/4 . M = imply dynamo runs. The Prandtl 
number Vr = u/rj for all runs was set equal to 1. The 
last columns give the measured parameters IZe = (u 2 }2 jvk u 
the Reynolds number and TZ\ = (u 2 ) 2 /u{(S7 x u) 2 ) 2 the 
Taylor Reynolds number. 1Z\ = (u 2 )5A/f (with A 2 = 
(u 2 )/((V x u) 2 )) the Taylor Reynolds number. The last col- 
umn gives the range of values obtained by fi = E M /E K the 
ratio of magnetic to kinetic energy. The simulations named 
A# were carried out in a grid of size 512 3 while the simu- 
lations named B# were carried out in a grid of size 1024 3 . 
The third 1024 3 -grid simulation started with the parameters 
of B3' and run up to t ~ 24 turn over times t u = l/fc„(u 2 )2 . 
afterwards it was continued wit the parameters B3. 



ergy ratio p, for all the examined runs. The data cover 
more than two decades of the parameter fi. Three dif- 
ferent behaviors can be observed. First, over the range 
/x (0.5< /x <20) the energy dissipation is independent of 
/i. This implies that e follows the Kolmogorov scaling 
e oc uf/£ even when the large scale magnetic energy is 
twenty times greater than the turbulent kinetic energy. 
The data include runs that vary from fully helical to non- 
helical, and strongly magnetically forced to dynamo runs 
and for 1Z\ ~ 100 to 1Z\ ~ 300. Thus this result seems 
to be very general and robust in this range. 

At n larger than 10 two new branches appear. The 
results of run A 12 that is fully helical and strongly mag- 
netically forced is marked by triangles in figurc[2] For this 
run both magnetic and kinetic energy is concentrated in 
the large scales building helical structures with very small 
turbulent fluctuations The large scale magnetic and ki- 
netic energy both increase with time keeping their ratio 
(i fixed while the small scale fluctuations and the dissipa- 
tion rates saturate with time. As a result the normalized 
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FIG. 1. Kinetic and magnetic energy evolution as a function 
of time for the three 1024 3 runs. The vertical dashed line 
indicates when the forcing parameters changed from the ones 
for B3' to the ones for B3. 

dissipation rate £/U 3 ms k u decreases with time resulting 
in the behavior seen in figEl The dynamics here are con- 
trolled by magnetic helicity condensates, and despite the 
large Reynolds number they do not follow a turbulent 
scaling. 

Run A13 marked by circles in figure [2] is strongly 
magnetically forced in order to achieve large values of 
fi within the time limitations imposed by the computa- 
tional costs. This run although in agreement with the 
turbulent scaling for fi < 20 it transitions to the scaling 
e oc jU 1 / 2 as indicated by the dashed line in the figure. 
This scaling can be understood if we consider that the 
main mechanism for cascading the injected energy is not 
the velocity shear S u oc U rms k u but rather the magnetic 
shear Sb oc B rms kb that shreds Alfven-wavepackets as 
they travel along chaotic magnetic field lines. The re- 
sulting scaling is e oc SbU 2 oc [7 3 fc(,/W 2 that is observed 
in figure [2] Thus the large scale field rather than sup- 
pressing the turbulence cascade it enhances it. 

None of the runs showed a weak turbulence scal- 
ing that would have implied according to [1] the scaling 
e oc U A /BLk u oc /i -1 / 2 . There are few possible inter- 
pretations for this result. First, just like the case of the 
uniform magnetic field, the absence of weak turbulence is 
can be explained by the entrapment of energy in modes 
that vary perpendicular to the local magnetic field so that 
lx <C &\\ that makes the nonlinear coupling strong. An- 
other possibility is the lack of uniformity of the magnetic 
energy density: regions exist in space with weak local 
magnetic field where eddies can be stretched with no re- 
sistance from magnetic tension. If the cascade in these 
regions dominates the scaling of turbulence becomes in- 
dependent of Bl- A final possibility is that magnetic in- 
stabilities reveal themselves altering the scaling (like run 
A13) before the magnetic field becomes strong enough 
for weak turbulence to manifest itself. 

The energy spectra E m , Ek for the high resolution runs 
B1,B2,B3 are shown in figure[3l For the Bl dynamo run 
the kinetic energy spectrum shows a power-law scaling 
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FIG. 2. Top panel: The energy dissipation rate e normalized 
by Urmsku as a function of /i — E M /E K for all runs. B ms = 
^J2E AI and XJ rmB = *J2E M . The dashed lines indicate the 
scaling e oc Ur ms k u and e oc Ur ms B r mskb oc /x 1 . 



slightly less steep than the /c~ 3 / 2 prediction. A clear 
change of slope is observed at k ~ 8 where E m (k) be- 
comes larger than Ek(k) is also observed. The magnetic 
energy spectrum does not show a power-law scaling. The 
absence of a power-law scaling might be linked to the 
absence of large scale magnetic structures in randomly 
forced dynamo runs. The B2, B3 runs despite the one 
order of magnitude difference in show similar results. 
The kinetic energy spectrum is well fitted by a k~ 3 / 2 
power-law while the magnetic energy spectrum is best 
fitted by a k~ 5 / 3 law. The fact that different exponent is 
measured for the velocity and for the magnetic field in- 
dicates that probably a clear inertial range has not been 
reached yet even at I024 3 grid sizes. Two reasons affect 
the measured spectrum exponents. First, at small wave 
numbers the large scale field Bl contaminates the iner- 
tial range making the spectrum look steeper. Second at 
large wave numbers the energy spectrum can be altered 
by bottle-neck effects that are far from well understood 
in MHD. The insets in figure [3] show the total energy flux 
U(q) = (u<(u ■ Vu - b ■ Vb) + b<(u • Vb - b • Vu)) 
(solid lines) and the energy flux due to the vortex stretch- 
ing term U u (q) = (u<(u- Vu)) (dashed line). u<, b< 
are the filtered velocity field containing only wavenum- 
bers with |k| < q. In all three cases the vortex stretching 
term plays a minor role in cascading the energy. 

The top panel of figure H] shows the measured power- 
law exponents from all the runs. The exponents were 
calculated by fitting power-law solutions Ek ~ A\ t k~^ K 
and E m ~ A m k~^ M in the range 2k u < k < where 
fc„ corresponds to the peak of the enstrophy spectrum 
k 2 Ek- The range of fitting can be seen in figGH by the 
vertical dotted lines. The bottom panel shows the dis- 
tribution of these exponents in the range 1 < n < 20, 
excluding run AI2 and the A1-A3 dynamo runs. Al- 
though the dispersion of these values is quite large the 
exponents do not show much dependence on amplitude of 
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FIG. 3. Kinetic and magnetic compensated energy spectra 
for the Runs B1,B2,B3 (from top to bottom). The straight 
lines indicates the scaling k~ 5 ^ 3 (solid) and fc~ 3 / 2 (dashed). 
The inset shows the energy flux IT(A;) (solid line) and H u (k) 
(dashed line). 



the magnetic field. Their values are concentrated around 
f or magnetic energy spectrum and fc~ 3 / 2 for 
the kinetic energy spectrum in agreement with the high 
resolution runs. 

Concluding, this work has showed an independence of 
the normalized energy dissipation rate (e oc U^ ms k u ) and 



spectral exponents (/3 M ~ 5/3 and /3 K ~ 3/2) on the am- 
plitude of the large scale magnetic held over a large range 
of fi and for a variety of forcing configurations. Nonethe- 
less, deviations were observed for very strong magnetic 
fields fi > 20 (for which (e oc B rms U^ ms kb) and for fully 
helical flows that formed magnetic helicity condensates 
in the large scales. Analysis of higher order moments Q 
and scale interactions [23[ might shed more light in the 
processes that control MHD turbulence. 
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